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Materials that exhibit loss or gain have a complex valued refractive index n. Nevertheless, when 
considering the propagation of optical pulses, using a complex n is generally inconvenient - hence 
the standard choice of real- valued refractive index, i.e. n 3 = Re(vr?)- However, an analysis of 
pulse propagation based on the second order wave equation shows that use of n 3 results in a wave 
vector different to that actually exhibited by the propagating pulse. In contrast, an alternative 
definition n c = ^/Re(n 2 ), always correctly provides the wave vector of the pulse. Although for 
small loss the difference between the two is negligible, in other cases it is significant; it follows that 
phase and group velocities are also altered. This result has implications for the description of pulse 
propagation in near resonant situations, such as those typical of metamaterials with negative (or 
otherwise exotic) refractive indices. 



I. INTRODUCTION 

Recent work in metamaterials and negative refractive 
index media 1 ll|El|lH,H,p,m,^] has focused attention 
on propagation in media with exotic values of permittiv- 
ity e and permeability /i, as well as those with significant 
loss or gain, where e and /i are complex valued. These 
material properties (i.e. e, /i) impact directly on the re- 
fractive index, and hence on the wave vector (3 and phase 
and group velocities ||, [T^] . 

When considering analytical solutions of the wave 
equation, it is often convenient to allow the propaga- 
tion wave vector (3 and refractive index n to be com- 
plex valued, based on the definition n 2 — c 2 efi, so that 
(3 = (uj 2 n 2 1 c 2 ) 1 / 2 . However, although this leads to many 
useful results, the approach also has some serious draw- 
backs. For example, the sign of the imaginary part of 
/3, which determines whether the wave experiences gain 
or loss, needs to be specified according to the chosen di- 
rection of propagation. Worse, in the envelope and car- 
rier description of pulse propagation, which is common in 
nonlinear optics (e.g. see |jll|), the presence of a complex 
wave vector in the carrier function is very inconvenient, 
since it requires the nonlinear coefficients to be adjusted 
to compensate for the distance propagated. In addition, 
determining other parameters such as the group velocity 
under these circumstances is also a non-trivial task (see 
e.g. iQ). For these and other reasons, it is often prefer- 
able to define a real-valued wave vector k and to treat 
the imaginary component separately. 

The standard approach is to simply define k as the 
real part of /3, i.e. k = ((jj/c)Me(\/n 2 ) = um s /c. 
However, an alternative definition based on k 2 — 
(w/c) 2 IRe (n 2 ) = oj 2 n 2 /c 2 has been used with advantage 
in studies of causality-based constraints for negative re- 
fraction jl3|, |l4| , although neither paper remarked on the 
non-standard definition. In that context, this alternative 



definition is required because it keeps the real and imag- 
inary parts of n 2 separate, and so ensures the Kramers- 
Kronig relations jl5| continue to hold, linking the two 
parts and enforcing causality. In contrast, the standard 
complex n is not required to be causal, although it is so 
in the case of passive (lossy) media (see e.g. [|l6[ p7[). 

In the present paper, the two definitions will be com- 
pared using the predictions of the second-order wave 
equation as the benchmark. It is shown that for field 
propagation in media with loss ( "passive" ) or gain ( "ac- 
tive"), where the use of a complex wave vector is particu- 
larly problematic, the alternative definition has the clear 
advantage that it exactly matches the spatial oscillations 
of the field. In contrast, the standard definition gives an 
imperfect match, and the description only recovers the 
true propagation due to the presence (and inconvenience) 
of additional correction terms. Note that the alternative 
definition (for n c ) is not in any sense equivalent to one 
based on an effective refractive index, such as might oc- 
cur in (e.g.) waveguides: it is an alternative choice of 
definition for the bulk refractive index. 

Because I focus on the propagation of waves, in sec- 
tion H I present a short description of the second order 
wave equation. Then, in section III, I give some defini- 
tions required for the handling of both the standard case 
(section IV) and the new alternative definition (section 
[V]) . After discussion of the similarities and difference be- 
tween the definitions in s ectio n VI, I end by presenting 



my conclusions in section VH 



II. THE SECOND ORDER WAVE EQUATION 

The second order wave equation is commonly used in 
optics (at least as a starting point) in descriptions of 
propagation, and results from the substitution of the V x 
H Maxwell's equation into the VxE one in the source- 
free case (see e.g. |ll[). In homogeneous media, with 



V 2 = di 



d 2 



d 2 and d a = d/da, the frequency space 



wave equation is 
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1 Also commonly called negative phase velocity (NPV) media 



V 2 E + (3 2 E = 0. 



(1) 



1 



2 



Here (3 2 = eficu 2 is the square of a complex propagation 
wave vector, since both e and fi can be complex. We can 
relate it to a complex refractive index squared quantity 
with 



0* 



(2) 



When considering the propagation of helds, it is conve- 
nient to split /3 2 up into two parts (e.g 
imaginary parts). Here I write f3 2 = k 2 
eqn. (Q) becomes 



its real and 

9 

^7 



V 2 E- 



k 2 E + ij 2 E = 0. 



(3) 



When considering this wave equation, we will usually 
want the first two terms to give plane- wave solutions, 
with the rest component containing loss and nonlinear- 
ity 2 . This is an important step, since although we might 
solve linear problems using a complex valued n, realistic 
situations are not so easily handled. 

The first two terms in eqn. (^), taken in isolation, 
have plane- wave solutions if k is real- valued; I call this 
the "underlying propagation" . The third term in eqn. (|^) 
is the "residual" component, which controls the discrep- 
ancy between the true propagation and the underlying 
propagation. Although in the case of small loss or gain 
the residual component will be only a weak perturbation, 
the theory presented here is valid for any strength. 

As an aside, if we specialize to the case of fields prop- 
agating along the z direction, using the carrier and en- 
velope models of pulse propagation |n], [H], ^l], p2| , 
we would write E(z,t) = A(z, t) exp[i(ujt — kz)] + c.c. 
to accommodate the rapidly oscillating behaviour of the 
carrier frequency: this carrier represents the underlying 
propagation for a specific frequency. This then leaves 
only the (usually) slowly varying envelope A{z, t), which 
would be affected only by the residual component. 

Returning to the wave equation of eqn. (|3|) , and tak- 
ing propagation along the z-axis, we can now factorize it 
using Greens functions Jll| [H], |23| , to give two first-order 
equations that are coupled only by the residual compo- 
nent. At the same time we can split the field into forward 
(E + ) and backward (£L) parts (i.e. set E = E + + E_), 
to give a pair of coupled, counter propagating, first order 
differential equations. These are 



d z E± = ±ikE± T ^ (£ 4 



E_) 



(4) 



Here the underlying propagation is, as desired, plane- 
wave like, since the first RHS term just adds an ikz be- 
haviour onto the frequency dependent iujt. The prop- 
agation is then modified by the second RHS term, i.e. 
the 7 2 -dependent residual component. A feature of this 



approach is that we see that any contribution (whether 
linear or not) that is included in the residual component 
will couple the forward and backward fields together (see 
M, |T) for more discussion). Since such terms are scaled 
by k in eqn. (^), they change (but in a simple way) under 
my alternative form for the refractive index. 

Here I consider only the one dimensional linear case, 
where 1 is independent of the field. This covers the cases 
of both loss and/or gain (i.e. in passive and/or active 
media); however for simplicity I will often only refer to 
loss; nevertheless the case of gain is always allowed for 
(since gain can be seen as "negative loss" ) . 

If we take the propagation to be of the form E + = 
Eq exp [i(uit — k'z)], with E- = 0, then eqn. (Q) gives 
us 



— ik! = —ik 



2fc' 



(•5) 



so that 7 2 < corresponds to loss for a forward propagat- 
ing wave. Further, if we consider instead the oppositely 
propagating wave, eqn.(Q) automatically ensures the nec- 
essary change of sign to ensure a loss stays loss, and a 
gain stays a gain. In contrast, when using a complex- 
valued n, care must be taken to ensure the correct sign 
(see e.g. M). 



III. DEFINITIONS 

We have that /3 2 and n 2 are (in general) complex val- 
ued, and lo and c are strictly real valued. Thus when 
choosing the propagation wave vector we need to decide 
what to do about the imaginary parts. Our choice then 
affects the performance, utility, and convenience of the 
refractive index, phase velocity, and group velocity. 

I now define some useful intermediate quantities to 
express the refractive index conveniently; I introduce 



and the angle 4> — Arg(n ) so that 



, 2 - r,?*'* 



n = rige 



4/2 



(6) 
(7) 



Whether or not specific values of cf> correspond to a neg- 
ative refractive index or negative phase velocity can be 
determined from the criteria for e and [i given in ]25[| 3 . I 
also define a reference wave vector k„ such that 



(8) 



The standard form for a real valued refractive index 



is 



no cos • 



(9) 



We can even incorporate diffraction in the rest by including the 



transverse parts of V 
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Note that the <j> used here corresponds to </>+ in the summary in 



3 



I have already noted that many treatments leave n as 
a complex valued quantity, leading to a complex wave 
vector k; and that while useful in many circumstances, in 
the context of pulse propagation it brings some significant 
disavantages. 

An alternative definition for the refractive index is 



e(n 2 ) = n Q ^/cos </>, 



(10) 



where n 2 satisfies the Kramers-Kronig relations [ |l5[ in 
partnership with the imaginary part Im (n 2 ) ; this defini- 
tion has already been used in the literature (e.g. see the 
recent H,[l4)). 



IV. THE STANDARD FORM 



If we rewrite eqn. ( [lq ) to incorporate the correction into 
the leading term, we get 



d z E± = ±ik s 



l-^sin^ 
2fc? 2 



E± T t-^- sin 2 | E_ 



E-). (17) 



As before, the first term on the RHS is gives plane-wave- 
like propagation, but now with a wave vector that differs 
from k s . 

I will now express the effective propagation wave vector 
in terms of k n and <f>. To simplify the description, I apply 
the usually excellent |^7j approximation that the effect of 
on the propagation can be ignored (i.e. set E- = 0). 
Hence, 



The standard form for the wave vector based on the 
standard form of refractive index (see eqn. (^])), 



fc = 



CO 



= ki cos 



kg — fan cos ■ 



(ii) 

(12) 



Thus k s is always real-valued, and can be negative in 
some circumstances. The phase velocity is then the usual 



c/n s , and the (inverse) group velocity simply v 1 



dk 
dui ' 



Let us now consider how this standard form of fc 2 looks 
when substituted into the second order wave equation. 
To do this let us express /3 2 in terms of fc 2 and fc 2 , 



[3 2 = k 2 s +zk 2 nl 2 s . 
with the residual behaviour described by 



'7.; 



sin (f> -j- % sin — 



(14) 



This standard choice of fc : 
wave equation of the form 



fc s leads to a second order 



d 2 z E + k 2 s E +ik 2 nl 2 s E = {). 



(15) 



When factorized, as briefly described in section |J, wc 
get a pair of coupled, counter-propagating, first order 
equations. These are 



d z E ± = ±ik s E ± T J-7s (E+ + E_) 



(16) 



Since the residual component tj 2 on the RHS of eqn. 
( |l6| ) contains a real part as well as an imaginary part, it 
is not pure loss. The real part will impose oscillations on 
the field as it propagates, thus altering the wave vector 
away from the assumed value fc s . However, the real part 
is quadratic in </>, being tx sin 2 ^, so for small losses the 
correction to the underlying propagation will be small. 



8 Z E+ = +ik' s E + - lysine/) E+, (18) 



with 



k s — k n cos — 



1 tan 2 - 



For <fi <C 1 , we then find that 

fc' c 2 ~ fc, 2 cos < 



(19) 



(20) 



Thus although I began with the standard definition, 
which assumes that the (forward-like) field will propa- 
gate with a wave vector fc = fc s , we see instead that it 
propagates with a wave vector fc ~ fc„^/cos 4>. As we will 
see, this approximation to the effective propagation wave 
vector is usually close to that of the alternative form dis- 
cussed below; the difference (for small loss) is of order 



(13) 4 



The standard phase velocity v p is 



k 2 



(21) 



However, if we were to use the effective propagation 
wave vector k' s we would get a different answer; in the 
case of the approximate form of eqn. (pfj|), it turns out 
the same as the alternate form given in the next section. 

The standard group velocity v a can be derived using 



Ikedukc — fc e 



— {d w n ) - {d u 4>) tan - + - 
uq 2 uj 



(22) 



Hence 



dujk s 



id 



"0 



- (d u 4>) tan - 



(23) 



Just as for phase velocity, if we were to use the effec- 
tive propagation wave vector k' s , we would get a different 
answer; in the case of the approximate form of eqn. (|2"c|), 
it turns out the same as the alternate form given in the 
next section. 
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V. THE ALTERNATIVE FORM 

The alternative form for the wave vector, based on the 
product e/i, (i.e. the square of the refractive index, see 
eqn. (Oil)), is 



fc 2 = — Re (n 2 ) = fc 2 cos < 



k„ yj COS (f). 



(24) 
(25) 



Thus k c is either real-valued or is pure imaginary. Real 
values of k c correspond to a regime of propagating waves, 
imaginary values to that of evanescent waves. The phase 
velocity is then u p — c/n c , and the (inverse) group veloc- 
ity simply Ug 1 = both will differ from the standard 
v p ,Vg, and are given below. Note that k 2 is related to k 2 s 
by 



k 2 



fc 2 cos <j> 
kl cos 2 t 



= 1 - tan J -. 



(26) 



With this alternative choice, it is simple to express (3 2 
in terms of our wave vector fc 2 , 



[3 2 = k 2 + k 2 n7 2 , 
with the residual behaviour described by 



*7 C 



!Sin< 



2 • 2 V 

l ls + sm 2" 



(27) 



(28) 



For small <p <C 1 , 7s and 7 C differ only by terms of order 
(j> 2 . Note that the loss-like part of the residual component 
(i.e. of Im(7 2 ) or Im(7 2 )) is the same for either form; 
but that only this alternative form of k (i.e. k c ) ensures 
that the residual component is purely lossy, and will not 
change the spatial oscillations of the field away from those 
of the propagation wave vector. However, the alternative 
form of k leads to the underlying propagation becoming 
evanescent if Me(n 2 ) < 0. 

With this choice of wave vector (i.e. k = k c ), the 
second order wave equation can be written 



d 2 E + k 2 r E 







(29) 



When factorized, as briefly described in section [n| we 



get 



fc 2 

8 Z E± = ±ik c E ± T 2^7 C sin (E+ + E_ ) . (30) 



The phase velocity u p is now faster than for the stan- 
dard definition, being 



or 



P f~2 



1 - tan 2 



(31) 



The corresponding group velocity u g can be derived 
using 



Ikrd/jjkr — fc„ 



no 



< tan i 



(32) 




FIG. 1: Comparison of k values, as a function of <j> = Arg(e^t). 
The causal choice k c is shown using a solid line when it is 
real valued, and dotted when imaginary ("ife c "); the standard 
choice (k 3 ) is given by the dashed line, with the approximate 
corrected form (k' a ) from eqn. (GOI) shown dot-dashed. 



Hence 



— dujk c — 



H (du,n ) 

n 



— (e? w 0) tan <j> 



(33) 



Here the comparison of u g with the standard form v g 
is less simple than f or phas e velocities: the prefactors dif- 
fer are different (^/ cos(0) compared to cos^); also the 
bracketed terms differ slightly (with tan0 not tan^). 

However, for <fi < ir/2, ^/cos(0) < cost), so that the 
group velocity u g is faster than the standard v g . 



VI. DISCUSSION 

As already noted, for small losses the standard and 
alternative definitions of n (and also those of fc) nearly 
coincide, but they diverge as the loss increases. Indeed, 
for (e.g. strongly resonant) situations where Ke(n 2 ) < 0, 
the underlying propagation (i.e. that defined by k s or 
k c ) can be of a completely different character. 

The simplest case is the trivial one where where 
Im(n 2 ) = 0. Here fc 2 = fc 2 , and both are always positive; 
both 7 2 and 7 2 are zero. The descriptions are identical. 

Next we add a small imaginary part to n 2 , with 
\4>\ <C 1, so that k s and fc c no longer match. The loss-like 
part of the residual component is (as always) the same in 
both cases, but a standard (fc s ) description will be modi- 
fied by an additional oscillation, giving an effective wave 
vector comparable to k c . This is perhaps the most typi- 
cal regime for device operation; being either the low loss 
case of normal (positive phase velocity) propagation, or 
the low loss case of NPV propagation. 

As <j) increases, the two descriptions diverge, as sum- 
marized on fig. 0. We see that the standard descrip- 
tion (fc = k s ) gives qualitatively similar behaviour for all 
1 01 < ?r; being one of a wave vector k s with added loss 
and a correction to achieve the true propagation wave 
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vector. Obviously, the larger the 4>, the larger the wave 
vector correction. 

The alternative choice of k = k c behaves differently. 
When \cf>\ = 7r/2, i.e. when n 2 = Im(n 2 ), the wave vec- 
tor k c vanishes, giving no underlying oscillatory evolu- 
tion as the field propagates. The only evolution is that 
given by the residual component, i.e. the loss speci- 
fied by Im(n 2 ). Then, as \c/>\ increases further, so that 
Me(n 2 ) = Re(c 2 e/i) < 0, we find that k c takes on an 
imaginary value: this is just the case of plasmons, where 
Re(e) G (—00, 0], but Rc(/i) 6 [0, 00). Here the imaginary 
k c means that underlying propagation becomes evanes- 
cent; and any loss then acts in addition to that. 

Note that the loss in the alternative description is sim- 
ply Im(n 2 ) - it differs from that used in the standard pic- 
ture. In particular note that this is not identical to the 
sum of the permittivity-based "loss" (i.e. Im(e)) and the 
permeability-based "loss" (i.e. Im(/i)). Further, at least 
in the case of doubly passive media |2(|, Im(n 2 ) < is in 
fact a criterion for NPV; i.e. loss is a criterion for NPV. 
More general statements on this relationship have been 
made when placing causality-based constraints on neg- 
ative refractive index media using the Kramers-Kronig 
relations g [L§. 

Lastly, whichever choice of k or n we make, it depends 
only on the sum of the complex phases of e and fi. In 
contrast, the summary given by J2(| shows that the NPV 
criteria of p5|] also depends on the difference of those 
phases. This sensitivity arises because the presence of 
NPV depends on the relative phases of the electric and 
magnetic fields; however the second order wave equation 
does not distinguish between the electric and magnetic 



responses, considering only their nett effect on the se- 
lected field (here, the electric field E). 



VII. CONCLUSION 

Here I have shown that the standard definition for a 
real valued refractive index (i.e. n = n s = Re(v / ?^ 2 )) 
is only an approximation to the true real valued refrac- 
tive index seen by a propagating optical pulse. Instead, 
the true propagation wave vector is based on the alter- 
nate definition n = n c = ^yWc{n 2 ). This conclusion was 
reached by examining how fields are actually propagated 
by the widely used electromagnetic second order wave 
equation, in the case where when loss (or gain) is treated 
as a modification to an underlying propagation based on 
a real- valued refractive index or wave vector. Treatments 
of pulse propagation that use this alternative n c (and 
hence k c ) will not only be using wave vector that ex- 
actly matches the propagation, but adjustments to that 
propagation will involve only gain or loss. In contrast, 
for the standard treatment based on n s , k s corrections to 
the spatial oscillation of the fields must be applied along 
with those for gain or loss. 
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